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Abstract 



HERA-B is a hadroproduction experiment located at DESY in Hamburg, 
Germany. The experiment produces B mesons and baryons by inserting thin 
wire targets into the halo of the proton beam circulating in the HERA stor- 
age ring. The B decays are studied to search for evidence of CP violation 
and constrain the angle (3 and possibly 7 and a of the CKM unitarity tri- 
angle. The experiment also produces B s mesons; these decays are studied 
to measure or constrain the mass difference Am s and width difference AT S 
between the two Bg/B s mass eigenstates. Finally, the large number of B's 
produced allows HERA-B to search for rare and forbidden decays such as 
B — > K^ifl^. The experiment is scheduled to begin running in early 2000. 
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1 Introduction 



HERA-B is a hadroproduction experiment that produces b mesons and 
baryons by colliding protons circulating in the HERA proton storage ring 
with fixed wire targets. The nominal interaction rate is 40 MHz. Down- 
stream of the targets are the following detectors: a silicon-strip vertex de- 
tector used to reconstruct decay vertices; a large-aperture (250 x 160 mrad 2 ) 
dipole magnet and thirteen drift chamber stations used to measure momen- 
tum; a Ring-Imaging Cerenkov counter (RICH) used to identify pions, kaons, 
and protons; an electromagnetic calorimeter used to measure electron and 
photon energies; and a set of gas ionization chambers interleaved with iron 
plates used to identify muons. More details about these detectors can be 
found in the HERA-B Technical Design Report 

Two main triggers are used for the experiment: a dilepton trigger and a 
high-p T trigger. The former identifies electron and muon pairs based on hits 
in the electromagnetic calorimeter or muon chambers. The hit positions are 
used as seeds for a Kalman-filter tracking algorithm which searches regions 
in wire tracking chambers for corresponding hits. When hits in several 
successive chambers are found, a track is reconstructed. For events with two 
or more tracks, the invariant mass of all opposite-sign track pairs (m i+i _) 
is calculated and events with m i+i _ above a threshold are passed. The 
threshold is chosen to efficiently select J ftp — > £ + £~ decays. The high-p T 
trigger works in the same manner except that hits in a set of pad chambers 
located within the magnet determine the track seeds used by the Kalman 
filter. The hit patterns that initiate a trigger correspond to tracks that 
have high transverse momentum (p T ) with respect to the beam. The mass 
threshold used for the high-p r trigger is chosen to efficiently select B° — > 
vr+vr-, B° -> K + n' , and 5° -> K + K~ decays]] 

Events that pass either the dilepton or high-p T trigger are input to a 
second-level trigger running in a "farm" of dedicated processors. This trigger 
projects the two dilepton or high-p r tracks into the silicon vertex detector, 
finds corresponding hits on silicon strips, and uses these hits to reconstruct 
silicon-based tracks. These silicon tracks are well- measured and are required 
to form a vertex some distance downstream of the nearest target wire. 

To measure CP violation, HERA-B measures the proper time distribu- 
tions dF(B° — > f)/dt and dT(B — > f)/dt and constructs the asymmetry 
parameter A CP (t) = (dT/dt — dT / dt) / (dT / dt + dT/dt). For the final state 

2 Charge-conjugate modes are included unless otherwise noted. 
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A CP {t) = -D T D M am(20) sm(xt) , (1) 

where D T is a "dilution" factor arising from mistagging the B° or B° flavor 
of the decay, D M is a dilution factor caused by tagging signatures from 
neutral B mesons which oscillate before they decay, (3 is the interior angle 
of the CKM unitarity triangle, x is the B°-B mixing parameter Am/T, 
and t is the decay time in units of B lifetime. Alternatively, HERA-B can 
count the total number of decays (corrected for acceptance) and construct 
the ratio of partial widths A CP = (T — T)/(T + T). This quantity equals 
—D T D M sin(2/3) £(t ), where £(t) = [sm(xt) + x cos(xt)]/(l + x 2 ) and t is 
the minimum value of decay time that a decay vertex must satisfy. Such a 
requirement is necessary to sufficiently reduce backgrounds. 

The proton interactions at HERA-B lead to a production asymmetry be- 
tween B° and B of a few percent; this must be measured and corrected for 
in order to detect a CP asymmetry. The correction is made by constructing 
the ratios 




where N B _ ¥n and N—^_ are the numbers of B° and B decays to a (copious) 
channel not exhibiting CP violation, e.g., B° — ► J/ifjK*° — ► t + £~ K + ir~ . 
The ratio (R-R)/(R + R) then equals (r - T)/(T + T). 
For N B ^ n » N B ^j, the uncertainty on sin(2/3) is: 

where N B is the number of B° — > J/tpKg candidates, 5 is the number of 
true signal events, B is the number of background events, and e tag is the 
tagging efficiency. The factor D T ^e^ indicates the performance of the 
tagging method. There are four methods to be used in HERA-B to tag the 
flavor of the decaying B° or B : 

1. the charge of a lepton with p (p T ) > 5 (0.8) GeV/c that originates 
from the semileptonic decay of the other b hadron in the event; 

2. the charge of a kaon with 5 < p < 50 GeV/c that originates from the 
b — > D — > i^r ± decay of the other 6 hadron in the event; 
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3. the "jet charge" defined by: 

Qjet = 9i I Pi • a| / S IP* " a l ' ( 3 ) 

i i 

where a is the direction vector of the B candidate as defined by the 
positions of the interaction and decay vertices, and i runs over all 
tracks in the decay vertex; 

4. the "same-side" tag, which is the charge of a track within a small rj-<p 
cone around the B° momentum vector which has minimum p T relative 
to the combined (B° + track) momentum vector. 

The first three methods have been studied extensively using a Monte Carlo 
simulation, and their expected efficiencies, dilution factors, and overall fac- 
tors D M D T ^/e^ are listed in Table |l[ 



Method 


e ta g D M D T D M D Ty/ /e tag 


Pi (Pt) > 5 (°- 8 ) GeV/c 
5 < p K < 50 GeV/c 

X)?ilpV a l / J2 IpV & I 

i i 


0.15 0.53 0.20 
0.48 0.29 0.20 
0.99 0.12 0.12 


All combined: 


0.29 



Table 1: The expected efficiencies, dilution factors D M D T , and overall fac- 
tors D M D T ^e^ for three tagging methods to be used in HERA-B. 

2 Measurement of sin(2/3) and sin(2o;) 

To estimate the precision with which HERA-B can measure sin(2/3) using 
B° — ► J/ipKg decays, we insert values for all parameters into Eq. (^). These 
values are listed in Table § the resultant uncertainty is Asin(2/3) = 0.17 
in one year of running. The precision with which HERA-B can measure 
sin(2a) using B° — ► tt + tt~ decays is also estimated using Eq. @; inserting 
values for the parameters (see Table |3|) gives A sin(2a) = 0.36 in three years 
of running. This estimate does not include theoretical uncertainty, which 
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is discussed below. The branching fraction for B — > ir + Tr~ is taken to be 
0.47 x 10~ 5 ||, and the signal-to-background ratio is taken to be one. Both 
estimates above include a lifetime cut on the decay vertex of 0.5r B , which 
is approximately 10 times the expected resolution in the vertex z position. 
The precision obtained for sin(2/3) would be a significant improvement over 
the current result from CDF||; the precision obtained for sin(2a) would 
constrain this parameter to a range much smaller than that allowed by 
current data [|[]. 





12 nb 




13 mb 


a bb/ ,(J inel 


9.2 • 10~ 7 


interaction rate 


40 MHz 


66 rate 


37 Hz 


fraction reconstructed as B° — ► J/t/jKg 


3.8 • 10~ 6 


(for / (66 -> B°X) ~ 0.7) 




reconstructed B° — > J/ipKg per 10 7 s 


1400 


statistical factor l/£(r B ) 


1.5 


(1/D M ) mixing of tagging 1? 


1.2 


£ + if + E ?i tagging factor (-D T ^e^) _1 


3.5 


A sin(2/3) after 10 7 s 


0.17 



Table 2: The expected uncertainty in sin(2/3) after one year of running. 

In addition to B° -> //^-Kg and B° — > 7r + 7r decays, there are other 
modes which provide sensitivity to sin(2/3) and sin(2a). These are summa- 
rized in Tables |3] and |B[ Each table has six columns: the left-most column 
lists the decay mode of interest; the second column lists the branching frac- 
tion for the decay as either measured from data or estimated in Ref. 
The third column lists the number of B + B decays reconstructed per year. 
This number is obtained by multiplying the branching fraction by the bb 
production rate, the fraction of 66 pairs which yield a B meson, the factor 
10 7 seconds of running per year, and the estimated HERA-B trigger and re- 
construction efficiencies. The fourth column lists the first-level trigger which 
accepts the decay. The fifth column lists the number of vertices arising from 
the decay chain and the number of kinematic constraints that can be used to 
reject background. The sixth column shows a sketch of the decay topology. 
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°bb 


12 nb 




13 mb 


a bb/ ^inel 


9.2 • 10~ 7 


interaction rate 


40 MHz 


bb rate 


37 Hz 


fraction reconstructed as B° — ► ir + ir~ 


5.6- 10~ 7 


(for f(bb -» ~ 0.7 and = 0.47 x 10~ 5 ) 




£?° — > 7r + 7r~ decays per 10 7 s 


210 


statistical factor 1/£(t b ) 


1.5 


(1/D M ) mixing of tagging B 


1.2 


£ + K + T,Qi tagging factor (Dy^/e^)- 1 


3.5 


A sin(2o;) after 3 X 10 7 s 


0.25 


(assuming S/B pa 1) 


(0.36) 



Table 3: The expected uncertainty in sin(2a) after three years of running. 
For simplicity, we assume the ratio of penguin amplitude to tree amplitude 
to be very small. 

These last two columns indicate the level of background suppression one 
expects for the decay mode: two or more decay vertices and a large number 
of kinematic constraints should provide substantial background suppression. 

Experimentally, the modes listed in Table |] are very promising: all have 
two decay vertices and at least five kinematic contraints. The branching frac- 
tions are substantial and indicate that HERA-B should reconstruct > 3000 
of these decays per year. Theoretically, the decays provide a clean mea- 
surement of sin(2/3): the dominant amplitudes have a single weak phase 
(Arg V* s V cb ) , and thus direct CP violation is negligible and Eq. (|T|) has neg- 
ligible theoretical uncertainty. 

The decay modes listed in Table || are more challenging: they have only 
one decay vertex and fewer kinematic contraints. Theoretically, the ampli- 
tudes may receive significant contributions from penguin diagrams which 
have a different weak phase from that of the tree diagram (ArgV^V^); in 
this case Eq. (||) also has a cos(xt) dependence and, in addition, the mea- 
sured value of a is shifted from the true value by an amount which depends 
on the unknown ratio of the penguin amplitude to the tree amplitude. This 
shift can be determined via an isospin analysis ||, but such an analysis re- 
quires measuring small branching fractions of difficult-to-reconstruct decay 
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modes such as B + — ► ir + n and -B — ► 7T 7T . Another method to de- 
termine a in the presence of penguins uses SU(3) flavor symmetry and the 
decays B° -> vr+vr", 5° ^+tt", and 5+ tt+K . Although the last de- 
cay will not be measured in HERA-B (it is not selected by the trigger), only 
the branching fraction is needed and this is expected to be well-measured 
at CLEO III|], BaBar§, §, and BELLE (u|. The three-body decays 
B° — > pir — > 7T + 7r _ 7r can also determine a in the presence of penguins via 
a Dalitz plot analysis [11]; however, the branching fraction for B° — > p°ir° 
is probably too low for this method to be useful in HERA-B. The modes 
B° — > pp which have two vector particles in the final state require an angular 
analysis to extract sin(2a), and this in turn requires a relatively large event 
sample. 

3 Measurement of sin(27) 

A number of decay modes have been discussed in the literature which provide 
sensitivity to the angle 7, which is the third angle of the CKM unitarity 
triangle. Some of these modes are listed in Table ||. The first three were 
proposed by Gronau and Wyler [^] ; by measuring their branching fractions 
and those of their charge conjugates, one can construct two triangles with a 
common base. The angle between the sides representing the B + — > D^ P K + 
and B~ — > D^ P K~ amplitudes is 27. Because only tree diagrams contribute 
to these decays, there is no uncertainty arising from penguin diagrams or 
final-state rescattering.P| Unfortunately, the expected branching fraction for 
B —> D K is too low for HERA-B to make a meaningful measurement. 



The next three modes were proposed by Dunietz [14] and are analogous to 
the first three but pertain to B° decays; their branching fractions are also 
expected to be too low for HERA-B. 

The amplitudes for B° — > K + n~ and B° ^ K~ir + can be used in con- 
junction with that for B + -> ir + K° (measured at CLEO / BaBar/ BELLE) 
to construct two more triangles and determine sin(27). This method^] 
requires estimating a color-allowed tree amplitude, and this introduces theo- 
retical uncertainty. However, the branching fraction for B° — > K + tt~ is large 
enough for HERA-B to collect a few hundred events per year, and such a 
sample could provide a first constraint on this angle. It is estimated [|l(| that 

3 However, if the D° is detected via a hadronic decay such as D° K~tv + , there is 
uncertainty arising from interference with a doubly-Cabibbo-suppressed amplitude such 
as D° -> K--K+ -see Ref. 111. 
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in one year of running, HERA-B could measure a difference in the rates for 
B° -> K + tt- and B° -> ^71"+ with a precision A.A CP /.A CP = 0.08 (0.16), 
for a signal-to-background ratio of unity (0.10). 

A bound 7 of the form < 7 < 7 or 7 < 7 < 180° which may have less 
theoretical uncertainty [17] can be obtained from the inclusive measurement 
B(B° — > K~ir + + B° — > K + ir~); however, this bound depends on the ratio 
B(B° -> if-7r+ + B° K+ir-)/B(B + -> K vr+ + B~ being 
less than one, which may not be the case (e.g., see Ref. Q). 

The decay B° ^ D* + tt~ is also promising, as the branching fraction is 
substantial and 2000-3000 of these decays should be reconstructed per year. 
Because the final state is not a CP eigenstate, fitting for the decay time 
dependence yields sin(2/3 + 7) rather than sin(27), which is advantageous if 
7 is near 90° (as possibly indicated by data from CLEO g |||). This decay 
proceeds via the b — ► cud transition; because all quark flavors in the final 
state are different, only tree diagrams contribute and there is no theoretical 
uncertainty arising from penguin diagrams or final-state rescattering. The 
ultimate success of this method and the others discussed above depends on 
as-yet-unknown backgrounds. 



4 B® Physics 

In addition to B d mesons, HERA-B also produces B s mesons; these allow the 
experiment to measure or constrain parameters of the B^-Bg system. For 
example, the decays B° -> Dj£+u £ , B° -► D~ir + , B® -> D s -vr+7r-vr+, B a s -► 
J/ijjK , and B s ° -» ir+if" can be used to measure the mass difference Am. 
or decay width difference AT S between the two B®/B s mass eigenstates. 
The mode B® — > DjK + (analogous to B° -> D* + ir~, i.e., no penguin 
contribution) can be used to measure 2(f) M + 7, where (f> M is the B®-B s 
mixing phase (expected to be very small). The mode B® — » J/ifxf> should 
not exhibit CP violation at HERA-B's level of sensitivity Jl9[ , and thus 
observing a CP asymmetry would indicate physics beyond the Standard 
Model. For all these modes, the number of events HERA-B expects to 
reconstruct per year is listed in Table ^. For most modes, the number 
of reconstructed events is > 100 per year, or > 300 events in three years. 
Unfortunately, the number of reconstructed B® — > D~K + decays is too low 
for HERA-B to constrain 2(p M + 7; this sample needs to be large in order 

to resolve the rapid B^-Bg oscillations. 
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The "reach" of HERA-B for measuring the mass difference Am s has 
been estimated based on the amplitude fit method of Ref. |2C| 1. This reach 
is taken to be that value of Am s in which the quantity S/N (signal to noise 
ratio) equals 1.645, where S/N is given by: 

| = ^1 fAl-2v)e-^ Am)a/2 . (4) 

In this expression, n is the number of B® or decays reconstructed, f s is 
the fraction which are true signal, rj is the probability of a mistag, and o t is 
the decay time resolution. Typical values for these parameters are listed in 
Table |8] and imply that in three years of running, Am s could be measured 
or constrained up to values > 24 ps _1 . This is superior to the current world 
limit Am s > 14.3 (95% C.L.) |lj. 

HERA-B can search for CP violation in the B® system by measuring 
the decay time distribution of B® — > K + K~ decays and looking for a de- 
viation from an exponential time dependence.^ If the time distribution is 
exponential, then it can be used along with that for B® — » D~tt + to de- 
termine the width difference Ar s [^]. Since dN KK /dt cx e _r i* = e~ r £*, 
and dN Ds7T /dt oc e -( r i+ r 2 )*/2 cosh(Ar s /2)t « e -(r 1 +r 2 )t/2 = e -{v L +v H )t/2 j 

T Ds7T - f KK = (T L + T H )/2 - T L = ATJ2. HERA-B expects to recon- 
struct ~ 35 Bg — > K + K~ decays per year, and a preliminary study indicates 
that such a sample could be sensitive to a difference Ar s /r s ~ 0.18. This 
sensitivity is superior to that of the current measurement from CDF 24]. 



5 Summary 

In summary, the physics potential of HERA-B is substantial. We expect the 
experiment to do well measuring sin(2/3), with an expected uncertainty of 
Asin(2/3) = 0.17 after one year of running. The experiment can constrain 
sin(2a), but this requires a better theoretical understanding of the penguin 
contribution to this mode. If the penguin contribution were very small, 
HERA-B could obtain Asin(2a) ~ 0.36 after three years of running for a 
signal-to-background ratio of unity. The experiment can potentially con- 
strain the angle 7 via B° — > K + ir~ and B — > D* + ir~ decays; the viability 

4 A non-exponential decay time distribution demonstrates interference between the 
B Sy H and B s ,l amplitudes to K + K~; thus there exists an interference term in the rate 
\{K+K-\H\B°)\ 2 which changes sign under CP. Such a term is CP violating - see e.g., 
Ref. @. 
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of these methods depends on unknown backgrounds. HERA-B can constrain 
Am s to be > 24 ps -1 in three years of running, which is an improvement over 
the current world limit of Am s > 14.3 (95% C.L.). HERA-B will also search 
for rare and forbidden dilepton decays such as B — > K*£ + £~, B® — ► (ft£ + £~, 
and B -» K^^e^ [||. Finally, HERA-B may find new phenomena such 
as measurable CP violation in 5° — > J/ipcj) decays. 
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Decay mode 



BR Rec. Evnts Trigger Vertices/ 
(xlO 5 ) per year Constraints 



Topology 



{J /iff -► 



7T + 7T 



(j/y> -► 

(if* if°7T°) 

(Y>' -► J/V'vr+vr-) 

(J/V> -»• £+0 



(■)// — > J/^ir + Tr~ 

(J/V> -► £ + r ) 
(if* - 



(Xcl ~> J 1^1) 



(j/y> -»• £+0 

(if* -► if^TT ) 



3.7 



1.9 



0.87 



0.58 



1.1 



0.52 



1400 



700 



330 



220 



430 



200 



fi/jb/ee 



/jL/jt/ee 



/jL/jt/ee 



2/5 



2/7 



2/6 



2/8 



2/6 



2/8 



Table 4: Methods to measure the CKM angle (3 in HERA-B. 
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Decay mode 



BR Rec. Evnts Trigger Vertices/ Topology 
(xlO 5 ) per year Constraints 



B° 
B° 



7T+7T" 



7T 



K+ 



0.47 
1.4 



210 
260 



high-p T 
high-p T 



1/2 
1/2 



B° — p+TT- 

B° — p°7r° 



4.4 
1.0 
0.1 



< 1900 


high-p T 


1/4 


< 420 


high-p T 


1/4 


< 42 


high-p T 


1/4 


<C 1200 


high-p T 


1/6 


< 27 


high-p T 


1/4 




B° 



p+p- 



B° -> p°p° 



2.9 
0.064 




Table 5: Methods to measure the CKM angle a in HERA-B. A "<" sign 
indicates that the trigger and reconstruction efficiencies assumed are opti- 
mistic. 
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Decay mode 



BR Rec. Evnts Trigger Vertices/ 
(xlO 5 ) per year Constraints 



Topology 



B~ -» D°K~ 

(D° -► K-TT+) 



B~ -» D"K~ 

(D° -► K+tT 



1.5 



0.0077 



0.23 



< 650 high-pp 



< 3 



high-pp 



< 98 high-pp 



1/4 



1/4 



1/4 



0.023 



< 10 high-p T 



2/5 



0.0077 



< 3 



D° cp K*° 



B° 

{D% P -► X+if 
(DOp^vr+vr-) 



0.0035 



< 1.5 



high-pp 



high-pp 



2/5 



2/5 



TV 



D*+tv- 

K-TT+) 



1.4 



7.4 



260 



high-pp 



nO . „0 t>-0 

B s^ P K S 



0.003 



< 3100 high-pp 



<C 1.5 high-pp 



1/2 

2/5 

2/5 
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Table 6: Methods to measure the CKM angle 7 in HERA-B. A "<" sign 
indicates that the trigger and reconstruction efficiencies assumed are opti- 
mistic. 



Decay mode 



BR Rec. Evnts Trigger Vertices/ 
(xlO 5 ) per year Constraints 



Topology 



B° s D-7V+ 

(D- -► (j)ir-,K* K- 

(4> -► 



12 



£> s 7T+7T 7T+ 



s° - D-e+u £ 



B° S ^J/^K 



*o 



r*0 



32 



320 



0.52 



< 280 high-p T 



< 740 high-p T 



900 high-p T 



100 



1/5 



2/5 



1/4 



1/4 



5° -» D~K+ 

s s 



0.32 



0.63 



< 7.4 high-p T 



< 15 high-p T 



1/5 



1/5 



g° -> K + K~ 



1.5 



35 



high-p r 



1/2 



(0 -► A"+A"-) 



5.5 



110 



\x\xjee 



1/4 



Table 7: 5^ decay modes to be collected in HERA-B. These modes can be 
used to measure Am s and Ar s , to constrain s'm(2(p M + 7), and to search 
for new physics. All D~ final states are reconstructed via D~ — > 4>ir~ — > 
K+K-7T- or Dj -► A"*°if- -► K+tt-R-. A "<" sign indicates that the 
trigger and reconstruction efficiencies assumed are optimistic. 
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D 


£> S -7T + 
s 7T + 7T 7T+ 


J/ijiK*" 


n (events/yr) 


~900 


~400 


~100 


f s (signal fraction) 


0.35 


0.6 


0.9 


^mistag (P rob - of mistag) 


0.24 


0.33 


0.33 


cr z (vertex resolution) 


600 fim 


450 /jm 


300 urn 


<J t (proper time resolution) 


0.11 ps 


0.081 ps 


0.054 ps 


((p B ) « 100 GeV/c) 








Am s reach in one year 








(S/N = 1.645) 


12 ps" 1 


13 ps" 1 


14 ps" 1 


Am s reach in two years 


14 


17 


21 


Am s reach in three years 


15 


18 


24 




(World limit: 


Am s > 14.3 ( 


3) 95% C.L.) 



Table 8: Decay modes to be used to measure or constrain the mass difference 
Am s , and the lower limits obtained in one, two, and three years of running. 
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